Relaxation of a Moving Contact Line and Landau-Levich Effect 
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The dynamics of the deformations of a moving contact line is formulated. It is shown that an advancing 
contact line relaxes more quickly as compared to the equilibium case, while for a receding contact line there is 
a corresponding slowing down. For a receding contact line on a heterogeneous solid surface, it is found that a 
roughening transition takes place which formally corresponds to the onset of leaving a Landau-Levich film. 



When a drop of liquid spreads on a solid suface, the contact 
line, which is the common borderline between the solid, the 
liquid, and the corresponding equilibrium vapor, undergoes a 
rather complex dynamical behavior. This dynamics is deter- 
mined by a subtle competition between the mutual interfacial 
energetics of the three phases, dissipation and hydrodynamic 
flows in the liquid, and the geometrical or chemical irregular- 
ities of the solid surface [|l[j. 

A most notable feature of contact lines, which is responsi- 
ble for their novel dynamics, is their anomalous elasticity as 
noticed by Joanny and de Gennes For length scales below 
the capillary length, which is usually of the order of 1 mm, 
a contact line deformation of wavevector k, denoted as h(k) 
in Fourier space, will distort the surface of the liquid over a 
distance |fc| . Assuming that the surface deforms instanta- 
neously in response to the contact line, the elastic energy cost 
for the deformation can be calculated from the surface tension 
energy stored in the distorted area, and is thus proportional 
to \k\, namely E cl = ^ J |§|fc||/i(fc)| 2 , in which 7 is the 
surface tension and 9 is the contact angle 

The anomalous elasticity leads to interesting equilibrium 
dynamics, corresponding to when the contact line is perturbed 
from its static position, as studied by de Gennes M\. Balanc- 
ing —^f- and the dissipation, which for small contact angles 
is dominated by the hydrodynamic dissipation in the liquid 
nearby the contact line, he finds that each deformation mode 
relaxes to equilibrium with a characteristic frequency (inverse 
decay time) u(k) = c\k\, in which c = -f9 3 /(3r)£) where 7/ 
is the viscosity of the liquid and I is a logarithmic factor of 
order unity ^ . The relaxation is thus characterized by a dy- 
namic exponent z, defined via ui(k) ~ \k\ z , which is equal 
to 1. The linear dispersion relation implies that a deforma- 
tion in the contact line will decay and propagate at a constant 
velocity c, as opposed to systems with normal line tension 
elasticity, where the decay and the propagation are governed 
by diffusion. This behavior has been observed, and the lin- 
ear dispersion relation has been precisely tested, in a recent 
experiment by Ondarcuhu and Veyssie [Q]. 

As an interesting example for nonequilibrium cases, corre- 
sponding to when there is an overall relative motion between 



the liquid and the solid, Landau and Levich studied the wet- 
ting of a plate vertically withdrawn from a completely wet- 
ting liquid at a constant velocity — v The case of partial 
wetting, where the liquid has a finite contact angle 9 e at equi- 
librium, has been studied by de Gennes [HQ]. He finds that 
the velocity v and the dynamic contact angle 9 are related as 
v = c(9^/6 — l)/2, and thus argues that a steady state is 
achieved in which the liquid will partially wet the plate with 
a nonvanishing dynamic contact angle for pull-out velocities 
less than c, while a macroscopic Landau-Levich liquid film, 
formally corresponding to a vanishing 9, will remain on the 
plate for higher velocities, as depicted in Fig. 1 [H], Note that 
at the transition there is a jump in the "order parameter" 9, 
from 9 e /\fi to zero. 

In the presence of defects and heterogeneities in the sub- 
strate, which could be due to (surface) roughness or chemical 
contamination, a contact line may become rough because it 
locally deforms so as to find the path with optimal pinning 
energy [Q]. This is in contrast to the case of a perfect solid 
surface, where the contact line is flat. The roughness can be 
characterized as a scaling law that relates the statistical width 
W of the contact line to its length L, via W ~ Z>, with the so- 
called roughness exponent £ being equal to 1/3 for the case of 
surface disorder with short-range correlations [^]. The contact 
line is also pinned by the defects, which means that a nonzero 
(critical) force is necessary to set the contact line to motion, 
through a depinning transition [Q] . It is also important to note 
that there may be numerous metastable states for the contact 
line due to the random disorder, leading to hysteresis in the 
contact angle [Bjjjj], 

Here we study the dynamics of the deformations of a 
nonequilibrium (moving) contact line. Using a balance be- 
tween the hydrodynamic dissipation in the deformed moving 
liquid wedge and the rate of change of the interfacial energies, 
we find that the dynamical relaxation is described by 

d t h(k,t) =-(c-v)\k\h(k,t) (1) 
in which v is the average velocity of the contact line, and 
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X(q,k-q) = -(2v-c)q(k-q)+3c\q\\k-q\ + (c-v)\k\(\q\ + 
| k — q\ — | k\) is the leading nonlinear correction. As compared 
to the equilibrium case (v = 0), we thus find that the linear re- 
laxation is faster for an advancing contact line (v < 0), while 
it is slower for a receding one (v > 0). In particular, for a 
contact line that is receding at the "terminal" velocity v = c, 
which coincides with the onset of leaving a Landau-Levich 
film, linear relaxation becomes infinitely slow and the domi- 
nant relaxation is thus governed by the nonlinear terms. 

We also take into account the effect of surface disorder on 
the moving contact line dynamics, which appears as a stochas- 
tic term in the right hand side of Eq. ([[]), and attempt to system- 
atically study the dynamical phase transition using renormal- 
ization group (RG) techniques. We find that the onset of leav- 
ing a Landau-Levich film formally corresponds to a roughen- 
ing transition of the (receding) contact line, which for a ran- 
dom substrate with strength g takes place at a critical velocity 
below c, corresponding to a dynamic contact angle 
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to the leading order. We combine our results with studies 
of the contact line depinning transition [H^Op, and propose 
a phase diagram for the system as depicted in Fig. 2. In par- 
ticular, we find that the phase boundaries corresponding to 
the dynamical phase transition and to the depinning transition 
meet at a triple point, and suggest that for stronger disorder a 
receding contact line will leave a Landau-Levich film imme- 
diately after depinning. 

Let us assume that the contact line is oriented along the 
x-axis, and is moving in the y-direction with the position de- 
scribed by y(x,t) — vt + h(x,t), as depicted in Fig. 3. If 
a line element of length dl = dxy/l + (d x h) 2 is displaced 
by Sy(x, t), the interfacial energy will be locally modified by 
two contributions: (i) the swept area in which liquid is re- 
placed by vapor times the difference between the solid-vapor 
75V and the solid-liquid 75^ interfacial energies, namely, 
{isv — lSL)dl5y/ y/T+ (d x h) 2 , and (ii) the work done by 
the surface tension force, whose direction is along the unit 
vector T that is parallel to the liquid-vapor interface at the 
contact and perpendicular to the contact line, as 7T • ydlSy. 
Note that we are interested in length scales below the capillary 
length, where gravity does not play a role. The overall change 
in the interfacial energy of the system can thus be written as 
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and we have made use 



in which a(x, t) = cos 1 • yj, 
of the Young equation: -fgy — jg^ = 7cos0 e . Note 
that both "forces" should be projected onto the y-axis when 
calculating the work done for a displacement in this direc- 
tion. For small contact angles and deformations one obtains 



a/I + (d x h) 2 cosa ~ 1 - 9(x, t) 2 /2, where 9(x, t) is the lo- 
cal contact angle. 

To calculate the dissipation, we assume that the contact 
angle is sufficiently small so that the dominant contribution 
comes from the viscous losses in the hydrodynamic flows of 
the liquid wedge [|lj] . For a slightly deformed contact line, we 
assume that the dissipation can be approximated by the sum 
of contributions from wedge-shaped slices with local contact 
angles 9(x,t). This is a reasonable approximation because 
most of the dissipation is taking place in the singular flows 
near the tip of the wedge [fflj^]. Using the result for the dis- 
sipation in a perfect wedge which is based on the lubrication 
approximation [p]JTT|], we can calculate the total dissipation as 



P = J dx^l + (d x h) 2 I 
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to calculate — ^ and P, and 
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We can now use Eqs.(Q) and 
set them equal to each other to derive the dynamical equation 
It yields 



2>T]i [v + d t h(x, t)] = 
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The above equation might simply be recovered by locally ap- 
plying the result of Ref. [^[ for straight contact lines, with the 
additional geometrical factor, which is needed when the direc- 
tion of motion is not perpendicular to the contact line, taken 
into account. 

To complete the calculation, we need to solve for the pro- 
file of the surface and the corresponding angles as a function 
of h. One can show that the surface profile z(x, y) near the 
contact line can be found as a solution of the Laplace equa- 
tion (d 2 + d 2 )z(x, y) = 0, so as to minimize the surface area. 
The solution reads z{x, y) = 9y+ J ^(3{k) exp (ikx — \k\y), 



where (3(k) = —6 



h(k) + f^\q\h(q)h(k-q)+0(h s ) 

is found using the boundary condition z(x, h(x)) = 
[§]. We therefore find 9{x) = 6{1 + J ^\k\h{k)e lkx + 
\j ^^f(k,k')h(k)h(k')e^ k+k '>] with /(A:, ft') = |ft + 
ft'|(lfc| + |fc'|) - (k + k 1 ) 2 + ftft', which can then be used in 
Eq.(Q) to yield Eq.(jl]), and the relation between v and 9 as 
described above and depicted in Fig. 1 . 

The linear relaxation of a moving contact line thus takes 
place at characteristic frequencies which obey the modified 
dispersion relation u(k) = (c — v)\k\, and as the onset of 
leaving a Landau-Levich film corresponding to v = c is ap- 
proached, the relaxaxtion becomes progressively slower. 

In addition to dissipation and elasticity, the dynamics of 
a contact line is also affected by the defects and hetero- 
geneities in the substrate, which are present in most practical 
cases. If the intefacial energies -fsv an d Isl are space depen- 
dent, a displacement of the contact line is going to lead to a 
change in energy as SE C \ = J dxg(x,vt + h(x,t))5y(x,t), 
where g(x,y) = jsv(x,y) - jSL(x,y) - (tsv - Tsl)- 
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Incorporating this contribution in the force balance leads to 
a noise term on the right hand side of Eq.(jl|) of the form 
r)(x,t) = j^g(x,vt) to the leading order. Note that this 
is a good approximation provided we are well away from 
the depinning transition, and the contact line is moving fast 
enough [p]J^,ppC|] . Assuming that the surface disorder has 
short range correlations (so that the correlation length is a mi- 
croscopic length a) with a Gaussian distribution described by 
(g(x, y)) = and (g(x, y)g(x', y')) = g 2 a 2 S(x - x')5(y - 
y'), we can deduce the distribution of the noise as: (r/(x,t)) = 
and ( V (x, t)rj{x',t')) = £j^5(x - x')S(t - t'). 

In the presence of the noise, the contact line undergoes dy- 
namical fluctuations. These fluctuations can best be charac- 
terized by the width of the contact line, which is defined as 
W 2 (L, t) = j- J dx(h(x, t) 2 ). Using the scaling form of the 
two-point correlation function, one can show that W ~ 
for intermediate times, while it saturates to W ~ Is at long 
times. Similarly, we can study the fluctuations in the order 
parameter field S8(x, t) = 6(x, t) — 9, and find 



(S6(x,t) 2 ) ~ l-B/t 2 ^/*, 



(6) 



where B is a constant. Note that the order parameter fluctua- 
tions approach a finite limit at long times provided f < 1. 

Keeping only the linear term in Eq.(|l]), we calculate the 
width of the contact line as 
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and, similarly, the order parameter fluctuations as 



(66(x,t) 2 ) ~l-B/t 2 
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for t ^> t — — — r . We thus find ( = and z = 1 within the linear 

(c— v) ^ 

theory. 

The nonlinear terms in Eq.(jj]) will modify the above results 
only when it becomes appreciable at long length scales, as 
compared to the linear term. The ratio of the two terms in 

Eq.(|l]) scales as ( c { ~!,)(L/a)< ~ L(t-v) ' and is thus appreciable 
only when the smallest time scale in the linear theory a/(c—v) 
becomes comparable to L/c. This happens near the onset of 
leaving a Landau-Levich film. 

Let us now attempt to systematically study the dynamical 
phase transition, corresponding to leaving a Landau-Levich 
film, using the RG scheme. The dynamical equation, which 
can be generally written in d dimensions as dth(k,t) = 

-u\k\h(k, t)-y ^X(q, k - q)h(q 7 t)h(k -q,t) + T){k, t), 

with X(q, k-q) = -Ai?(fc-g)+A2|?||fc-g|+As|fc|(|?|+|fc- 
q\ — \k\), belongs to the general class of Kardar-Parisi-Zhang 
equations [O]. We take a noise spectrum given by t)) = 
and (rj(k,t)r)(k',t')) = 2D{2Tr) d 6 d {k + k')5(t - t'), and 
employ standard RG techniques following Ref. to cal- 
culate the RG equations describing the flow of the coupling 
constants. We find 
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to the one-loop order, in which Sd is the surface area of a unit 
sphere in d dimensions. 

To study the fixed point structure of the above 
set of flow equations, it is convenient to intro- 
duce the dimensionless coupling constant U = 
5 d (7r/a) d+1 (Ai + A 2 )(A 2 + A 3 )L>/[2(27r) d ^ 3 ], and thus 
have z = 1 + U and £ = 1 — U at the fixed points. 
The flow equation for U reads: dU/dl = — (d + 1)U + 
[6 + (Ai + A 2 )/(A 2 + A 3 )] U 2 /2, which has two stable fixed 
points at U — (linear theory) and U — oo (strong cou- 
pling), as well as an intermediate unstable fixed point at 

j j _ j j* = 2(d+l) 

U U - 6+(A 1 +A 2 )/(A 2 +A 3 )' 

For U < U*, the nonlinearity is irrelevant and the ex- 
ponents are given by the linear theory, i.e. ( — and 
2 = 1, while for U > U* the behavior of the system is 
governed by a strong coupling fixed point which cannot be 
studied perturbatively. The fixed point at U* corresponds to 
a roughening transition of the moving contact line. The ex- 



ponents at the transition are z = 1 



2(d+l) 



and 
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2(d+l) 



e+{x 1 +x 2 )/(X2+X3) 
which are nonuniversal. The 



6+(Ai+A 2 )/(A 2 + A 3 )' 

strong coupling fixed point should presumably describe the 
Landau-Levich film. 

We can also study how the transition is approached by 
linearizing the flow equation near the fixed point. Setting 
U = U* + SU, we find dSU/dl = (d + 1)SU that would 
imply divergence of the correlation length near the transition 
as £ ~ |(5[/|~" with v = l/(d+ 1). The correlation length 
corresponds to the typical size of rough segments in the con- 
tact line, which should diverge as the transition is approached. 

Using the equation for the phase boundary U = U* with 
d = 1, and the values of the coupling constants that we 
can read off from Eq. ([[]), we map out the phase diagram of 
the system as depicted in Fig. 2. In particular, we find the 
asymptotic form of the phase boundary for weak disorder as 
reported in Eq.(0) above, and the nonuniversal exponents z 

17 , 15 /Hit •a/3/' 9 \2/3 nnr i f _ 5 15 ( Hit \ 1/3/ 9 

to the leading order, and v = 1/2. The order parameter fluc- 
tuations at the transition are given by Eq.(^) with the above 
choice for z and (, and remain finite since ( < 1. We also 
find that the triple point is located at 9 t /0 e = 0.887 and 
gt/(j6 2 ) = 0.138, which is interestingly still within the weak 
disorder limit. 

It is important to note that since our approach is only valid 
when the system is away from the depinning transition, it is 
not clear that we can trust the prediction of our theory in the 
vicinity of the triple point. In particular, how the two phase 
boundaries merge at the triple point, and the nature of the third 
phase boundary (dashed line) remains to be clarified. Further- 
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more, the above calculation is only restricted to the case of 
small contact angles where the dissipation is dominated by the 
singular hydrodynamic flow near the tip of the liquid wedge. 
Should other mechanisms of dissipation become important, 
the above picture would be altered [ pj[ ] . 

We finally mention that there could be two types of experi- 
ments to check these results, which we hope to motivate. The 
first type would be an analogue of the Ondarcuhu-Veyssie ex- 
periment [Q], which could probe the relaxation of a moving 
contact line and measure the velocity dependence of the dis- 
persion relation. 

The second type of experiments would correspond to a sys- 
tematic study of the onset of leaving a Landau-Levich film for 
receding contact lines on a disordered substrate. In particular, 
it would be interesting to look for a roughening of the contact 
line before the Landau-Levich film is formed. 
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FIG. 1. Dynamic contact angle as a function of pull-out 
velocity 




FIG. 2. The suggested phase diagram of a contact line on 
a disordered substrate. The depinning transition line is taken 
from Ref. [jiOQ, corresponding to the receding and the advanc- 
ing contact angles. The asymptotic form for the roughening 
transition line is given in Eq.(g). 




FIG. 3. The schematics of the system. 
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